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Optimal scale-free network with a minimum scaling of transport efficiency for 
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Average trapping time (ATT) is central in the trapping problem since it is a key indicator characterizing the 
efficiency of the problem. Previous research has provided the scaling of a lower bound of the ATT for random 
walks in general networks with a deep trap. However, it is still not well understood in which networks this 
minimal scaling can be reached. Particularly, explicit quantitative results for ATT in such networks, even in 
a specific network, are lacking, in spite that such networks shed light on the design for optimal networks with 
the highest trapping efficiency. In this paper, we study the trapping problem taking place on a hierarchical 
scale-free network with a perfect trap. We focus on four representative cases with the immobile trap located 
at the root, a peripheral node, a neighbor of the root with a single connectivity, and a farthest node from 
the root, respectively. For all the four cases, we obtain the closed-form formulas for the ATT, as well as its 
leading scalings. We show that for all the four cases of trapping problems, the dominating scalings of ATT 
can reach the predicted minimum scalings. This work deepens the understanding of behavior of trapping in 
scale-free networks, and is helpful for designing networks with the most efficient transport process. 

PACS numbers: 05.40.Fb, 89.75.Hc, 05.60. Cd, 89.75.Da 



I. INTRODUCTION 

As a paradigmatic random walk, trapping problem has 
received increasing interest within the scientific commu- 
nity, which was first introduced in the pioneering work 
by Montroll more than forty years ago 1 -. It is a kind of 
random walk in graphs in the presence of a single deep 
trap positioned at a given location, absorbing all particles 
(walkers) that visit it. Trapping process is relevant in a 
large variety of other dynamical processes occurring in a 
number of different complex systems. Frequently cited 
examples include light harvesting in antenna systems 2 ^—, 
energy or exciton transport in polymer systems^—, page 
search or access in the World Wide Web&iS, to name a 
few. Because of its practical relevance, it is of signifi- 
cant importance to address trapping problem in diverse 
complex systems. 

A fundamental quantity pertaining to trapping prob- 
lem is the trapping time (TT), often referred to as the 
mean first-passage time (MFPT)ii _ — . The TT for a 
node i is the expected time taken by a particle starting 
off from i to first visit the trap. The average trapping 
time (ATT) is defined as the average of TT over all pos- 
sible starting nodes in the system, which can be used as 
an indicator of the efficiency of trapping. In the past few 
years, trapping problem in different kinds of graphs has 
been intensively studied, including the square and cubic 
lattice s 16 ' 17 , dendrimersi£~— , the Sierpinski gaske t 21 ' 22 
and Sierpinski towe r 23 ' 24 , the T— fractal 2 ^— and its ex- 
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tensio n 28 ' 29 , as well as scale-free graphs 3 ^—. These stud- 
ies showed how the ATT scales with the size of systems 
with various topological properties, unveiling the non- 
trivial effects of graph structure on the behavior of the 
ATT. 

More recently, it has been proven 3 - 6 - that for trapping 
problem in any graph with a trap fixed at an arbitrary 
node, the possible minimal scaling for the ATT is propor- 
tional to the graph size and the inverse degree of the trap 
node, which is universal and provides a maximal scaling 
for the lower bound of ATT for trapping in an arbitrary 
network with an immobile trap. Thus, a network is called 
optimal if this minimal scaling can be achieved for any 
node as a trap. However, since the MFPT from one node 
to another depends on the source-target distance 3 - 7 -, for 
trapping in a graph, the fact that the minimal scaling can 
be reached when the trap is fixed on one trap does not 
necessarily lead to the conclusion that the minimal scal- 
ing can be obtained when the trap is located at another 
node. For example, for trapping in Cayley trees^ - as a 
model of polymer networks 3 -^—, the minimal scaling can 
be reached when the trap is at the core node, but it can- 
not be achieved if the trap is at a boundary node. Thus, 
it is interesting to design or find optimal graphs, where 
the minimal scaling for ATT can be reached for any trap 
node. In particular, it is useful to derive closed-form so- 
lutions to ATT for separate trap nodes having minimal 
scaling since they are instrumental to quantitative un- 
derstanding of theoretical models and to approximate or 
numerical solutions. 

In this paper, we study the trapping problem in a hi- 
erarchical scale-free network^ 3 -. We study four cases of 
trapping problems with the immobile trap being fixed at 
four representative nodes, i.e., the root node, a periph- 



eral node, a neighboring node of the root with a single 
degree, and a farthest node from the root, respectively. 
For all these four cases, we derive analytically the explicit 
expressions for the ATT, based on which we further ob- 
tain their leading scalings for large network sizes. We 
show that the four trapping processes are all very effi- 
cient since their ATT grows linearly with the network 
size for the worst case. Moreover, we show that for the 
four representative cases of trapping problems, the pos- 
sible minimum scalings for ATT can be achieved. In this 
sense, the network being studied is optimal for trapping 
processes, which is helpful for designing networks with 
the highest trapping efficiency. 



II. CONSTRUCTION AND PROPERTIES OF THE 
HIERARCHICAL SCALE-FREE NETWORK 

We first introduce the model of the hierarchical scale- 
free network^, which is constructed in an iterative man- 
ner. Let G g (g > 0) denote the network after g iterations. 
Initially (g = 0), the network Go contains one node with- 
out any edge, which is called main hub (or root) node. 
At iteration g = 1, to generate G\ we introduce two more 
nodes and link them to the original node in Go . The two 
new nodes are named peripheral nodes of Gi, and the 
root node of Go is also the root of G\. For g > 1, G g is 
obtained from G s _i by adding two new copies, denoted 

by G g _ x and G g 2 } 1 , of G 9 _i to the primal G g _i, with 
all peripheral nodes of the two replicas being linked to 
the main hub of the original G s _i unit . The hub of the 
original G g _i and the peripheral nodes of the two du- 
plicates of Gg-i form the main hub node and peripheral 
nodes of G 9 , respectively. Repeating indefinitely the two 
steps of replication and connection leads to the hierarchi- 
cal scale- free network. Figure Q] illustrates schematically 
the network G3. 

It is easy to obtain that the number of nodes in G s , de- 
noted by N g , is N g = 3 9 . According to the node degree, 
all these nodes can be separated into four different sets 
by applying the method in^i: the peripheral node set P, 
the locally peripheral node set P c (1 < c < g), the root 
node set H comprising only the main hub of G g , and the 
local hub set H c (1 < c < g). For detailed explanation, 
we refer to reference^. All nodes in a set have the same 
degree, and the degree for a node in sets P, P c , H, and 
M r is, 



and 



K P (g) = g, 
Kp, c (g) = c, 

K h (g) = 2(2" -1), 

i=i 



K h 4g)=Y,2 t = 2(2 C -1). 
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FIG. 1. (Color online) Iterative construction process of the 
hierarchical scale-free network. The red nodes are peripheral 
nodes. 



respectively. In addition, the number of nodes in each of 
these four sets is, respectively, 



2 C 3 9 ' 



and 



2 x 3 9 ~ 



(5) 



(6) 



(7) 



(8) 



Based on the above quantities, it is easy to derive that 
the average degree of all nodes is 



K(g) = 4 



2 V 
3 



(9) 



which is approximately equal to a constant 4 when g is 
large enough, implying that the network is sparse. 

In G g , the maximal distance from the root to any other 
node equals to g. Let ¥ g denote the set of those nodes 
in Gg at a distance g from the main hub. We call such 
nodes the farthest nodes of G g . As shown above, G g 
is composed of a primal G 9 _i and two replicas of G g _i, 

i.e., G g % and G g %. 



The farthest nodes of G g must be in 



the two subgraphs Ggj^ and Gj^. Specifically, ¥ g con- 
tains the farthest nodes of the primal central subgraph, 

i.e., Gg-2, forming G 



(2) 



(i; 



1 and G g 



(2) 



Then, the number 



3 



of nodes in ¥ g , denoted by |F S |, satisfies the following 
recursive relation: 



\¥ a \ = 2\F g _ 2 \ 

Considering the initial conditions | 
Eq. (fTU)) is solved to yield 



(10) 

2 and |F 2 | = 2, 



|F 9 | = 



2(s+ 1 )/ 2 

2 9/2 



g is odd, 
g is even. 



(11) 



The network being studied presents some typical prop- 
erties as observed in many real-life system o 43 ' 45 . It is 
scale-free 4 ^ with its degree distribution P(k) following 
a power law form P{k) ~ fc 7 , where the exponent 
7 = 1 + . It displays the small- world behavior— with 
its average distance increasing logarithmically with its 
size 4 ^. Moreover, the network has an obvious hierarchi- 
cal structure that has also been observed in real-world 
networks, e.g., metabolic networks 4 ^. It is the precur- 
sor, probably the first model for hierarchical scale-free 
networks. It is thus of great interest to study how the 
unique structure affects dynamical processes taking place 
on the network. 



III. TRAPPING IN THE HIERARCHICAL SCALE-FREE 
NETWORK 

After introducing the model and properties of the hi- 
erarchical scale-free network, in this section we study 
analytically a particular random walk — the trapping 
problem — occurring on G g , where a deep trap is posi- 
tioned at a certain node. 

We first consider the simple random-walk model on 
a general network with N nodes. At each time step, 
the walker jumps from its current position to any of its 
neighboring nodes with an identical probability. Assume 
that node j is trap node. Let Tij denote the trapping 
time for node i, i.e., MFPT from node i to node j, which 
is the expected time for a walker staring off from node 
i to first visit the trap j. The highly desirable quantity 
related to the trapping problem is the ATT, Tj, which is 
the average of Tj over all the N source nodes distributed 
uniformly over the whole network. By definition, Tj is 
given by 



T, 



1 N 

-Yt. 



(12) 



In the sequel, in order to explore the impacts of the 
trap position and its degree, we will study analytically 
four cases of trapping problem performed on the hierar- 
chical scale- free network G g , with the perfect trap located 
at the root node, a peripheral node, a neighboring node 
of root node having a single degree, and a farthest node, 
respectively For these four representative trapping prob- 
lems, we will explicitly determine the ATT and show how 
their dominating behaviors scale with the network size. 



A. Trapping with the trap being positioned at the root 
node 

We first address the case when the trap is fixed at the 
root node of G g . To this end, we define a multiple trap 
problem with all the peripheral nodes being the traps. 
Let Th(g) denote the ATT to the root node; and let the 
T p (g) denote the ATT when all peripheral nodes are oc- 
cupied by traps. 

In order to evaluate Th(g) and T p {g), we introduce two 
intermediate variables T p ^(g) and Th, p (g). The former 
represents the MFPT from an arbitrary peripheral node 
to the root node of G g , while the latter Th tP {g) stands 
for the MFPT from the root node to any of the 2 9 arbi- 
trary peripheral nodes in G g . In Appendix YK[ we derive 
analytically the two intermediate variables T p ^(g) and 
Th,p{g), which read 



and 



(13) 



(14) 



T Kv(9) = Ul) "I, 



respectively. It should be mentioned that Eqs. (fTBj) 
and (|14[) were first found ir* 3 ^ by using the approach of 
generating functions. The obtained results for T p ^{g) 
and Th tP {g) arc very useful quantities for the following 
calculations. 

Having obtained the intermediate quantities, we now 
determine the quantities Th{g) and T p (g). From the net- 
work structure, the following recursive relations can be 
established: 

Thig) = ln(g - 1) + \ [T p (g - 1) + T Pth (g)} (15) 

and 

T P (g) = I [T h [,g - 1) + T h , P {9)\ + \T p {g - 1). (16) 
After some algebra, Eqs. (fl~5j) and (JTS]) can be recast as 
3T h (g) - T h (g - 1) = 2 [T p (g - 1) + T Pih (g)] (17) 

and 

3T p (g) - 2T p (g - 1) - T h , p (g) = T h (g - 1). (18) 
From Eq. (|18[) . we can further have 

3T p (g + 1) - 2T p (g) - T hiP (g + 1) = T h (g), (19) 
which, together with Eq. (|T5|) . gives 

3 [3T p (g + 1) - 2T p (g) - T h , p (g + 1)] 
- [3T p {g) ~ 2T p (g - 1) - T h . p {g)\ 
= 3T h (g) -T h (g-l)=2[T p (g-l) + T Pth (g)] , (20) 



4 



namely, 



T p (g + 1) = T p (g) + - [3T htP (g + 1) - T h , p {g) + 2T p . h {g)] . 

(21) 

Inserting Eqs. (fT3|) and (fT4|) into Eq. ([21]) and considering 
the initial condition T p (l) = 1/3, we can solve Eq. ([2"Tj) 
to obtain the explicit expression for the ATT when all 
peripheral nodes are occupied by traps: 



T P (g) 



20 



9 



9 



19 



(22) 



which is consistent with the previous result in££. Substi- 
tuting Eq. (|2"2"|) into Eq. (fTT)) and considering the initial 
value T/i(l) = 2/3, we can solve Eq. (|T7|) to obtain the 
exact formula for the ATT to the root node: 



32 /3V 8 34 
Thi9) = Y\2 -9 9 ~-9- 



(23) 



Notice that Eq. ([2U[) agrees with the result derived in^ 2 . 
by using the approach of generating functions. 

We proceed to express T p (g) and Th(g) in terms of 
network size N g , with an aim to obtain the dependence 
of the two quantities on N g . Recalling N g = 3 9 , we have 
g = lniV 3 /ln3 and 2 9 = (iV a ) ln2 / ln3 , which allow to 
represent Eqs. (f22[) and (|23|) as a function of N g as 



and 



9 v ^ 9 i n 3 9 v > 



Th (g) = ^f N y-in2/ln3 

yy ' 9 V 9 ' 9 In 3 9 



* N > 34 (25) 



When the system is very large, i.e., N g — > oo, 

T h {g) ~ (JV s ) x - lll2 / ln3 (26) 

and 



( jy jl-ln2/ln3 i 



(27) 



both of which scale sublinearly with the network size. We 

note that Eqs. (|26j) and ([27)) have been perviously derived 
in 50 

by using another technique, which is different from 
the one adopted here. 

Since the degree of the root node is Kh(g) = 2(2 9 — 1), 
which is can be expressed in term of network size N g as 
K h (g) ~ (A s ) ln2 / ln3 when N g is very largei^, then T h (g) 
can be expressed as a function of N g and Kh(g) as 



T h {g) 



N a 



K h (g) 



(28) 



which grows proportionally to the network size and in- 
verse degree of the root node as the trap. 



B. Trapping with the trap being located at a peripheral 
node 



We now consider the second case that the trap is fixed 
at one of the peripheral nodes. As shown above, there are 
|P| = 2 9 peripheral nodes in G g . We label sequentially 
these peripheral nodes by 1, 2, • • ■ , 2 9 from left to right. 
For convenience, we classify all the peripheral nodes into 
g + 1 sets denoted by (0 < i < g): For i = 0, f3 — {1}; 
while for 1 < i < g, (3, = {x^' 1 < x t < 2 1 }. In addi- 
tion, let Bi(Q < i < g) be the union of the sets (3k with 
< k < i, namely Bi = \Jl=oPk- Since for this partic- 
ular trapping problem, the selection of trap position has 
no effect on the ATT, with loss of generality, we choose 
the node belonging to (3 as the trap and denote Tg (g) 
as the ATT. 

Prior to deriving the ATT Tjs (g), we define and de- 
termine some new quantities. Let Tp i+1 ^ i {g) denote the 
MFPT for a walker starting from an arbitrary node in 
/3j_l_i to an arbitrary node belonging to Bi in H g . Then, 
we have the following relation: 

Tp i+1 , Bi (9) - -E [l + T h , p (k)+T 0i+1 , Bi (g)] + 

i 9 r i fc_1 



k=i+l L 



l=i 



(29) 



The two terms on the rhs of Eq. ([29[) can be elaborated 
as follows: the first term explains the process that a 
walker starting from a node in takes one step to 
reach a local hub node that has no links to other pe- 
ripheral nodes except those in (3i+i, and then jumps 
Th,p{k) +Tp i+1) Bi{g) more steps to hit a target. The sec- 
ond term describes the process that the walker starting 
from takes one step to reach a local hub that has a 
link connected to peripheral nodes not in ft+i, then takes 
Th,p(k) + \ Xw=i Tft + i,8fe(<?) steps to visit a destination. 



From Eq. ([29[) . we can derive 



and 

r /Jt ,8i-i(s) = 



9-1 



k=i+l i=i+l 



(30) 



g-i + 1 



g ^3 k-1 



(31) 



which, together with Eq. ([T4^). leads to 

Tp^.M = ^ZT^ T ^M9)- (32) 

Using the initial condition Tp gt B g _ 1 (g) = 8(3/2) 9 — 8, we 
can solve Eq. ([3"2"j) to obtain 



Tp i+u Bi (g) 



4(3 9 - 23) 

{g-i)X 



(33) 
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The above obtained quantity Tp i+lt B i (g) enables us to 
obtain the ATT Tg Q {g). By construction, we have the 
following relation 



TbM= 3 



9-1 



T h (g - 1) + T Kp (g) + 1 £ 7> fe+lA (.g) 
fe=o 

1 9-2 

T p (g - 1) + T^ 9 ,b 9 _i (ff) + g E T ^ +1 >8 fc (f) 

fc=0 

I^^k^-ij+r^) + 1 g T , fc+1 , Bfc ( ff ) 

i=2 L fe=0 

1 

+T p (i - 1) + T ( g i , Bj _ 1 (ff) + ^^Tf) k+1 ^ k {g) 



(34) 



The three terms on the rhs of Eq. (|34|) can be accounted 
for as follows. Note that G g consists of three copies of 



H, 



g _x: central G g _i, G^ 1 ^, and Ggj^. The first term 
on the rhs of Eq. describes the contribution to ATT 
made by nodes in the central G ff _i. The second term 

presents the trapping time of nodes in the subgraphs 

(2) 

Gr_i- The third term is a little complicated, which ex- 



,(2) 



,(i) 



plains the trapping time of nodes in the subgraph G g _ 1 . 

Inserting the above obtained related quantities into 
Eq. we can obtain the explicit expression for the 

ATT T Bo (g) as 



T B (ff) 



1 

39 



16 



6 
9 

(fj 



3 

62 
~9 



9-1 1 

Y- 

i=l 6 
3X9 



70 
~9 



3x9 



j=i+l ■> 



6(3 9 - 29) 
29-H 



(35) 



It is not difficult to find that the term with the highest 
exponent occurs when i = 1 and j = g. Moreover, in the 
infinite network size limit, i.e., N g —> 00, we have 



T Bo {g)~ N g /\og 3 N gi 



(36) 



which, together with Eqs. (|27|) . indicates that when the 
trap is fixed at a peripheral node the trapping process 
displays a lower efficiency than the case when all periph- 
eral nodes arc traps. Thus, the number of traps sensi- 
tively affects the behavior of trapping processes taking 
place on the hierarchical scale-free network. Based on 
N g = 3 9 and Eq. (JTJ) , we have K p (g) = g = log 3 A 9 , 
which together with Eqs. (|36|) leads to 



(37) 



T Bo (g)~N g /K p (g), 



implying that T Bo (g) is proportional to network size and 
the inverse degree of the trap node. 



C. Trapping with the trap being placed at a root's 
neighbor with a single degree 

Here we address the third trapping problem with the 
trap being positioned at a neighbor of the root, which has 
only one edge. Note that among all neighbors of the root, 
there are only two neighbors having a single connectivity, 
denoted by x and y, respectively. For convenience, let x 
be the trap node. The quantity we are concerned with is 
the ATT to trap node x for trapping in G gi denoted by 
T x (g). In order to determine T x (g), we first determine 
the MFPT, Th, x (g), from the main hub to node x in G g 
By construction, Th tX (g) can be written recursively as 

T h ,M = 27^T) + 2(2^1) [1 + 1 + ThA9)] + 
1 9 

2( 29-i) E 2 ''[ 1 + T ^^ z )+ r ^(g)]- ( 38 ) 



The first term on the rhs of Eq. (|3"5]) describes the fact 
that the walker, starting from the main hub, requires only 
one time step to hit the trap node x. The second term 
explains the process by which the walker first jumps to 
node y in one time step, then takes one step to return 
to the root, and continues to jump Th yX {g) more steps 
to reach the target. The last term accounts for the fact 
that the walker first makes a jump to a peripheral node 
or local peripheral node belonging to P or (2 < i < g), 
then takes T p ^h(i) time steps to the hub, and proceeds to 
the node x, taking Th, x (g) more time steps. 

Equation (|38[) can be simplified to 



T h , x {g) = 1 +£2* + J22 1 T p , h (i). 



(39) 



i=l 



Considering the initial condition Th, x (l) = 3 and 
Eq. ([12]), we can solve Eq. (|39p to produce the exact 
solution for Th tX (g), which reads 



T hia; ( 5 )=4(3 9 -29)-l 



(40) 



Note that the expression in Eq. (j4H)) is actually a special 
case of the more general formula derived in a different ap- 
proach ir>22. Then, the quantity T x (g) can be accurately 
evaluated as 

1 N — 1 
T x (g) =T h {g)- — + — ^ — T hiX (g) 

(41) 

which can be expressed in terms of network size N g as 

32 



T x (g) = AN g - 4 (Ag ln2 / ln3 + ^(JV s ) 1 -" l2 / ln3 + 

9 



AfN Ni„ 2 /in3-i _ 81n7Vg _ 79 
[ 9 > 9 In 3 9 ■ 



(42) 
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Thus, when N„ 



oo, 



read 



T x (g) ~ N g , 



(43) 



that is, the leading term of T x {g) grows linearly with 
the network size but inversely proportional to the trap's 
degree, which is exactly equal to 1 in this particular case. 



D. Trapping with the trap being fixed at a farthest node 



Now we focus on the case when the trap is located 
at a farthest node in G g . Since in this case, the ATT 
does not depend on the location of the trap, without 
loss of generality, we choose the leftmost farthest node 
as the target that belongs to G g X \ and we let Tf(g) be 
the ATT for this special trapping problem. Below, we 
first concentrate on the MFPT from the main hub to the 
trap in G g , denoted by Thj(g), based on which we will 
determine Tf(g). Furthermore we will show that both 
Thj(g) and Tf(g) have the same leading scaling. 



Pg(i) 



>i+3 



- 1 



and 



hg{i) = 2 



9 V 2 



16 /3 



9—* 



9 



(44) 



(45) 



respectively. 

After obtaining the expressions of related quantities, 
we next determine the MFPT Thj(g) from the main hub 
to the leftmost farthest node, as well as the ATT Tf(g). 



2. Exact solution and leading scaling for the MFPT from 
the root to the leftmost farthest node 



In order to find the explicit formulae for Thj(g). We 
distinguish two cases: (i) g is odd and (ii) g is even. 

When g is odd, the target belongs to V\. In this case, 
we have 



1. Related definitions and quantities 

In order to derive the expressions for Thj(g) and Tf(g), 
we introduce some more variables. For those nodes of 
G g that belong to G g _ x or G g _ x , we can classify them 
in the following way. Let , H g -i(0 < i < g — 1) be the 
set of the local hub nodes that are directly connected to 
g — i classes of local peripheral nodes in P c , and let V g -i 
(0 < i < g — 1) denote the set of the local peripheral 
nodes whose neighbors arc g — i different local hub nodes 
belonging to H c . In addition, we assume that T-L g = H 
and V g = P. 

The specific structure of the network shows that for a 
particle starting from the main hub to one of the |F g | far- 
thest nodes, it must follow the path H g -^V g ^f Hg-i — > 



T, 



9-2 



H 



H 



g-(i+2) 



9-3 
-> • ■ 



' -> TV(t-l) -> Hg- 

T-L\ or V\. For the particular case 



Vg-(i+l) 



that the leftmost farthest node is chosen as the target, 
the path should be definitely as follows: each time the 
walker starting from a current local hub in Hg-i, its next 
goal must be a local peripheral node in Vg—u+x) belong- 
ing to the subgraph G g \ i+2 s that is one of the three 
components of G g -i, then it continues to jump to the 
main hub of a subgraph G g -^ + 2) that is in the central 
component of G g -i. In this way, the walker moves on 
until it reaches the trap. 

According to the above definitions, it is necessary to in- 
troduce two more quantities p g (i) and h g (i). The former 
is the MFPT from a node in V g -i to any of its neighbor- 
ing nodes in % g _u + -\\ and the latter is the MFPT from 
a node in H g -i to any of its neighbors simultaneously 
belonging to "P g „( i+1 ) and G g _, i+2 y In Appendix IB1 we 
provide the detailed derivation for p g (i) and h g (i), which 



T hJ (g) 



ThAa) + 2 T /3 a ,B s -i(f) 



k=0 



k=0 



By plugging Eqs. ([11 ]) . ([55 )l . and (|35]) into Eq. flU), 
we obtain a closed- form solution to T^j[g) given by 

W 9) = f<3»-,>-f(l)% 29+ ™ <«> 

When g is even, the target is in Hi. In this case, 
Thj{g) can be calculated by 



T h j(g) 



j2p a m + J2h g (2k+i) 



(48) 



k=0 



k=0 



Substituting Eqs. (flg j) . ([53 ]) . (|33 ]> . and (jlB ) into Eq. @gj|, 
after some algebra, Eq. (|48p is solved to yield the explicit 
expression for Thj(g), given by 



r„ (9) = |(3»-2»)-f 



+ 2g. (49) 



Equations (|47p and ([49]) can be expressed, respectively, 



in terms of the network size N g as 
2\nN 179 



In 3 + 15 



(50) 
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and 



Then, the ATT Tf(g) can be determined as 



Thj{9) = y [N„ - (N g ) h ^^] + ™ [{Ng) ^/^ i] 



21niV g 
In 3 



(51) 



both of which indicate that for large networks, i.e., N g —> 

oo, 

T hJ (g)~N g , (52) 
behaving as a linear function of the network size. 

3. Closed-form formula and dominating scaling for ATT 

Before evaluating the quantity Tf(g), we introduce two 
new quantities T p j(g, n) and Thj(g, n), which denote the 
MFPT to the leftmost farthest node for a walker starting, 
separately, from a local peripheral node in V g - n and a 
local hub node in 7i g - n , which is part of the walking path 
from the root node to the leftmost farthest node. These 
two intermediary quantities and the interesting quantity 
Tf(g) can be determined by distinguish odd and even g. 

When g is odd, we have 

.g-i _ i g — 1 _ i 

2 1 2 1 

T v j{9,n)= Pa( 2l )+ E M2* + l)> (53) 



since for odd g, only when n is even, the quantity 
T p j(g,n) is meaningful. Plugging Eqs. (|44)) . and (|45|) 
into Eq. (|53|) leads to 



o 5 



(54) 



Analogously, for Thj(g,n) we have 

3-1 _ 1 .9-1 _ I 

2 2 

Thj(g,n)= PaW+ E M2* + l), ( 55 ) 

* = 2" ^~ 2 

which, combining Eqs. (HU, and ([^5]) into Eq. gives 



^»)^(3»--)-f-(i) a -?(r s "" 

32 „„ 69 
+ y 2« + 2,-2n + -. 



r / ( fl ) = -[T ft ( ff -i)+T h>/ ( fl )] 

+\[T P {g-l)+T Pth {g) + T KS {g)] 



+ 37i E 2 xr ^ 2 [Tp(ff-2i-2) + r Pl/ ( ff ,2i)] 



+ £ [3 9 - 21 - 3 (T h ( 5 - 2i - 3) +T hJ (g,2i + l)) 



+3 g - 2l - s (T p (g -2i-3) + T p , h (g - 2i - 2) 



+T fc , / (fl,2t + l))] + 



(57) 



The three terms on the rhs of Eq. (JSHJ) account for, 
respectively, the contribution of nodes in primal G 9 _i, 



G g _i, and G y g l 1 that form G g . Using above obtained 
related quantities, Eq. ([53)) is solved to obtain 

656 /3\ 9 28 {2\ B , 10 , 907 

45" ' 

(58) 

which can be expressed in terms of network size N g as 

20 f, T ,, T a n 2/ln3l 656 ,„ r N l-l n 2/ln3 



-,(2) 



2>(ff) = f [iV 9 -(iY 9 ) ln2 / lna ]-^f (iV,) 1 - 



28 1+ln2/ln3 10 In TV 907 

+ ¥lnT + l5"- (59) 



Similar to the case of odd g, for even g we can obtain 

m , s 16 „«s 656 /3V 32 /2\ 9 10 934 

T / ( ff ) = y(3 9 -2 9 )-— y +y5+T ^ 

(60) 

and 



Tf( 9 ) = f [N g - (iV 9 ) ln2/ln3 J - ^ (/v,) 1 —/" 1 



656 



■?<".) 



-l+ln2/ In 3 



+ 



101n7V„ 934 



+ IF • (61) 



5 v " a/ ' 9 ln3 45 

Both Eqs. (|59j) and (|6Tj) show that for large networks, 



T f (g) ~ N g , 



(62) 



behaving as a linear function of the network size, a scaling 
identical to that of Thj(g). Notice that for odd and even 
g, the degree of the farthest node is 1 and 2, respectively. 
Thus, Tf(g) scales proportionally to the network size and 
the reciprocal of the degree of the farthest as the trap. 

The phenomenon that the leading term of Tf(g) dis- 
plays the same behavior as that of Thj(g) can be ex- 
plained based on the following heuristic arguments. Note 
that G g contains three subgraphs, each of which is a 
copy of G g -i. For those nodes in the central sub- 
graph, their MFPT to leftmost farthest node is equal 
to Th(g — 1) + Thj(g), the dominating term of which is 
(56) T hJ (g); for nodes in the subgraphs G^l 1: their MFPT to 
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the leftmost farthest equals T p (g — 1) +T p j l (g) +Thj(g), 
the leading term of which is also T^ j(^); while for nodes 

in the subgraphs G < ^} 1 encompassing the leftmost far- 
thest node, their MFPT to the target is less than Thj(g). 
Therefore, for all nodes in G g , the dominating term of the 
ATT Tf(g) is proportional to network size N„, which is 
similar to that of Tw(<?). 



E. Result analysis 

In the preceding text we have studied four represen- 
tative cases of trapping problems with the trap located 
at the main hub, a peripheral node, a neighbor node of 
the root with a single degree, and a farthest node from 
the root, respectively. We have shown that their trap- 
ping efficiency exhibits rich behavior. For the four cases 
of trapping problems, the leading term of ATT exhibits 
evidently different dependence on the network size. It 
can grow sublinearly or linearly with the network size, 
or behaves as a linear function of network size by a log- 
arithmic correction, which shows that the degree of trap 
node plays an important role in the trapping efficiency. 
Our results also demonstrate that for the four cases of 
trapping problems, the transport processes are very effi- 
cient with the ATT increases at most linearly with the 
network size. 

Although for the four trapping problems, the trapping 
efficiency displays distinct scalings with the network size, 
Eqs. (|28)l. (j37) . (|43]l , and ((62]) indicate that the leading 
scaling of ATT for all four cases grows inversely propor- 
tional to the degree of the trap, regardless of its position. 
For example, for the two cases of trapping problems when 
the trap is fixed on a peripheral node or a farthest node, 
the dominating scaling of ATT is identical. In fact, ex- 
tensive numerical simulations also verify that for all trap- 
ping problems in the hierarchical scale-free network with 
a single trap, as long as the degree of the trap is identical, 
the leading behavior for their trapping efficiency is also 
the same. 

It has been proved^ that for trapping problem in a 
general sparse network having N nodes with node j be- 
ing the trap, the scaling of the lower bound for ATT Tj 
varies with the network size N as Tj ~ N/dj, where dj 
is the degree of the trap node j. For the four cases of 
trapping problems in the hierarchical scale-free network, 
this minimal scaling for ATT can all be achieved. In 
most of previously studied networks, this minimum scal- 
ing cannot be reached. For example, for trapping in the 
(1, 3)— flower and the (2, 2)— flower with the same degree 
sequence, when the trap is located on a largest node, 
their ATT display distinct behaviors, but both are very 
larger than the minimum scalin g 51 ' 52 . In this sense, the 
hierarchical scale-free network being studied exhibits the 
most efficient configuration (the optimal structure) for 
random walks with a perfect trap fixed at a given node. 



IV. CONCLUSIONS 

Previous works have shown that for isotropic random 
walks in general networks in the presence of a single deep 
trap, the leading scaling for the least ATT is proportional 
to the size of the network and the inverse degree of the 
trap. In this paper, we have presented an in-depth anal- 
ysis on four particular cases of trapping problems in a 
hierarchical scale-free network, with the trap being lo- 
cated at the main hub, a peripheral node, a neighbor 
node of the main hub with a single degree, and a farthest 
node, respectively. For all these four cases, we have de- 
rived closed-form formulae for the ATT, as well as their 
dominating scalings that are all equal to the predicted 
minimal scalings. In this context, the network under con- 
sideration has an optimal structure that is advantageous 
to efficient trapping. 

Our work may have practical implications for designing 
networks, especially scale-free networks, where minimiz- 
ing the transport efficiency is a central goal. For exam- 
ple, with the emergency of new preparation techniques, 
now it is possible to prepare new synthesized polymeric 
materials with a very complex geometry, e.g., hierarchi- 
cal scale-free topology, which is favorable to transporta- 
tion and diffusive dynamics and thus can be utilized as 
potential artificial antenna systems for light harvesting. 
Moreover, our technique for computing the ATT is rela- 
tively general, which also applies to other networks. For 
instance, for trapping in dendrimers with a trap fixed on 
a peripheral node, the ATT can be easily determined by 
using our method. 
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Appendix A: Determination of quantities T p , h (g) and 

Th,p(a) 

According to the particular construction of the hierar- 
chical scale-free network, we can find that the two quan- 
tities T p _h(g) and Th tP (g) obey the following recursive re- 
lations: 



T P,h{g) 



9-1 



(Al) 



and 



T h, P {9) 



Ef=i % 



s-i 



2« + 5^2 i (i + r Pih (i) + r hJI ( ff )) 



; = 1 



(A2) 

The two terms on the right-hand side (rhs) of Eq. (|A1[) 
can be explained as follows. The first term is based on 
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the fact that the walker takes only one time step to first 
reach the root. The second term describes the fact that 
the walker first makes a jump to a local hub node be- 
longing to Hi, then takes Th lP {i) time steps, starting off 
from the local hub, to reach any node in P, and con- 
tinues to jump T p ,h(g) more steps to reach the root for 
the first time. Analogously, the two terms on the rhs of 
Eq. (|A2|) arc based on the following two processes. The 
first term describes the fact that the walker, starting from 
the root, requires only one time step to hit a peripheral 
node. The second term explains such a process that the 
walker, starting off from the root, first jumps to a local 
peripheral node belonging to Pi in one time step, then 
makes T p ^(i) jumps to the root, and proceeds to any 
node in P, taking Th, P (g) more time steps. 
Equations (|A1[) and (|A2[) can be recast as 



9-1 

Tp,h(d) = 9 + ^T Kp {i) 

i=l 



9-1 



(A3) 



Th, P {9) = 2 - -^r + ^ E 2 %,h®> ( A4 ) 



i=l 



and 



respectively. From Eqs. (|A3|) and (|A4[) , we obtain 



T Pth (g + l)=g + l + J2T hiP (i) (A5) 



and 



T htP (g + 1) = 2-1 + ^ T ^ 2%, h {i). (A6) 

i=l 

Equation (|A5[) minus Eq. (|A3[) leads to 

T p>h (g + 1) - T p . h (g) = 1 + T, hp (g). (A7) 
Similarly, Eq. (|A6|) minus Eq. (|A4[) times 1/2 yields 

Th, P (9 + 1) - \ThM = 1 + \T Pth (g). (A8) 

Applying the initial conditions T Pi ^(l) = 1 and 7/^(1) 
1 . we can solve Eqs. (| A7|) and (|A8[) to arrive at the ex- 
plicit formulas for T p _h(g) and Th, p {g), which are given 
by Eqs. (fT3"]) and HU, respectively. 

Appendix B: Computation of quantities p g (i) and h g (i) 

According to the structure of H g , p g {i) and h g (i) sat- 
isfy the following relations: 



Pg(i) 



1 



l+[hg(i-l)+Pg(i) + l] 



9~\ 

g-i-2 \ 

+ J2 [l + T h . p (k)+p g (t)} (Bl) 



and 



h 3 (i 



2(29 



J- - J 2»- < [l +PS - 1) + + i x 2«— 1 



[H-T^fc-i-l) +&,(*)] 



9-i-2 > 

+ Y, 2 fe [l + T p , h (fc) + ^«] . 
fe=i J 



(B2) 



Equation (|B1|) can be explained as follows. The first 
term describes the process that the walker starting off 
from a node in V g -i may directly go to the target. Al- 
ternatively, the walker can jump to the neighboring node 
in - H ff _(i_ 1 ), then takes time h g (i — 1) to reach one of the 
local peripheral nodes in V g -i, and continues to bounce 
p g (i) steps to hit the target; this process is explained by 
the second term. The last term represents the fact that 
the particle goes to a local hub in FJ C (1 < c < g — i — 2), 
from which it takes an average steps Th, p {k) to return to 
one of the local peripheral nodes in V g -i, and then moves 
on average p g (i) steps to arrive at the target. Analo- 
gously, we can explain Eq. (|B2I) . 

Merging similar terms of Eqs. (|B1I) and lB2[) leads to 



g-i-2 



p g (i) =h g (i-l)+g-i+ ^ T htP (k) 



(B3) 



k=l 



and 



1 9-1 

h g (i) = Ap g {i - 1) + 2" + T pA9 - * - 1 



2« 

g-i-2 



k=l 



29 



-i-2 



(B4) 



fc=i 



Substituting Eq. (|B4[) into Eq. (|B3|) . we can obtain the 
recursive relation for p g (i) as 

p fl (i) = 4p a (i-2) + — j — ^ 2 k +T Pjh (g-i) + 



on 

g—i—1 



k=l 



g-i-2 



2s 



^ 2 fc r M (fc) + (. 9 -*) + "^ r hJ ,(fc) 



k=l 



k=l 



(B5) 

Considering the initial conditions p g (0) = 64/9 x (3/2) 9 — 
7 and applying Eqs. ([13]) and $n§, Eq. (|B5)l can be solved 
to yield 



Pg(i) 



■)i+3 



3X9 



1 



9 V 2 



1. 



(B6) 



Plugging Eq. (|B6[) into Eq. (|B4[) , we obtain the analytical 
expression of h g (i) as 



_ oi+4 



16 f'd 
~9~ 



+ 3. 



(B7) 



fe=i 
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